SOME EXPERIMENTS WITH INTEGRAL APOLLONIAN CIRCLE PACKINGS 



ELENA FUCHS AND KATHERINE SANDEN 

ABSTRACT. Bounded Apollonian circle packings (ACP's) are constructed by repeatedly inscribing circles into the tri- 
angular interstices of a configuration of four mutually tangent circles, one of which is internally tangent to the other 
three. If the original four circles have integer curvature, all of the circles in the packing will have integer curvature as 
well. In | SI], Sarnak proves that there are infinitely many circles of prime curvature and infinitely many pairs of tangent 
circles of prime curvature in a primitiv^Jintegral ACP. In this paper, we give a heuristic backed up by numerical data for 
the number of circles of prime curvature less than x, and the number of "kissing primes," or pairs of circles of prime cur- 
vature less than x in a primitive integral ACP. We also provide experimental evidence towards a local to global principle 
for the curvatures in a primitive integral ACPs. 



1. Introduction 

Start with four mutually tangent circles, one of them internally tangent to the other three as in Fig.[T] One can 
inscribe into each of the curvilinear triangles in this picture a unique circle (the uniqueness follows from an old 
theorem of Apollonius of Perga circa 200 BC). If one continues inscribing the circles in this way the resulting 
picture is called an Apollonian circle packing (ACP). A key aspect of studying such packings is to consider the 
radii of the circles which come up in a given ACP. However, since these radii become small very quickly, it is more 
convenient to study the curvatures of the circles, or the reciprocals of the radii. Studied in this way, ACP's possess 
a beautiful number theoretic property that all of the circles in an ACP have integer curvature if the initial four have 
integer curvature. The number theory associated with these integral ACP's has been investigated extensively in 
HGLMWYL iFOl . and llKOl . 




Figure 1 . Packing Circles 



Key words and phrases. Number theory, computational number theory, diophantine equations. 
K. Sanden was supported by NSF Grant 0758299. 



1 



2 



ELENA FUCHS AND KATHERINE SANDEN 



A central theorem to any of the results in these papers is Descartes' theorem, which says that the curvatures 
(vi, v'2, V3, V4) of any four mutually tangent circles satisfy what is called the Descartes equation, 

2 2 2 2 2 

(1.1) F(vi,V2,V3,V4) = 2(vj + v 2 + v 3 + v 4 ) — (vi + v? + V3 + V4) = 0. 

where a circle which is internally tangent to the other three is defined to have negative curvature (see [Cx| for a 
proof). Given this formula we may assign to every set of 4 mutually tangent circles in an integral packing P a 
vector v e Z 4 of the circles' curvatures. We use Descartes' equation to express any integral ACP as an orbit of a 
subgroup of the orthogonal group Of(Z) acting on v. This subgroup, called the Apollonian group, is specified in 
|GLMWY|, and we denote it by A. It is a group on the four generators 
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derived by fixing three of the coordinates of v and solving F(\) = for the fourth. Note that each Si is of order 
2 and determinant —1. Also, Sj fixes all but the ith coordinate of v € Z 4 , producing a new curvature in the ith 
coordinate. 

In their paper [GLMWY | the five authors Graham, Lagarias, Mallows, Wilks, and Yan ask several fundamental 
questions about the curvatures in a given integer ACP, which have mostly been resolved in [F0|, |F1 1, |BF|, and 
|KO|. In particular, they make some observations about the congruence classes of curvatures which occur in 
any given ACP, and suggest a "strong density" conjecture, that every sufficiently large integer satisfying these 
congruence conditions should appear as a curvature in the packing. In MBFL the first author and Bourgain prove 
a weaker conjecture of Graham et.al. of this flavor that the integers appearing as curvatures in a given ACP make 
up a positive fraction of N. Proving the "strong density" conjecture would be significantly more difficult. In this 
paper, we use the /?-adic description of the Apollonian orbit from BF11 to formulate this conjecture precisely and 
provide strong experimental evidence in Section[3]in support of it. Our conjecture is specified further in the case 
of two different ACP's in Section|3] It is stated generally here. 

Conjecture 1.1. Local to Global Principle for ACP's: Let P be an integral ACP and let P24 be the set of residue 
classes mod 24 of curvatures in P. Then there exists XpGZ such that any integer x > Xp whose residue mod 24 
lies in P24 is in fact a curvature of a circle in P. 

We note that X above depends on the packing P under consideration. In this paper, we investigate two ACP's 
which we call the Bugeye packing Pp and the Coins packing Pq- These packings are represented by right action of 
the Apollonian group on (—1,2,2,3) and (—11,21,24,28), respectively (see Fig.|2]for a picture). In the case of Pp, 
our data suggests that Xp B exists and is < 10 6 , as we find no integers x > 10 6 which violate the above conjecture. 
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The data for Pq, however, suggests that Xp c exists, but that it is > 10 8 . Namely, there are integers x > 10 8 in certain 
residue classes in the set S24 which do not appear as curvatures in the packing we consider. We explain this further 
in Section[3] 




Another interesting problem regarding ACP's is counting circles of prime curvature in a given packing. Sarnak 
proves in jSjJ that there are infinitely many circles of prime curvature in any packing. In light of this, we give a 
heuristic in Section|2]for the weighted prime count \j/p(x): 

(1.3) Yp{x)= £ log(a(C)) 

a(C)<x 
a(C) prime 

where C is a circle in the packing P and a(C) is its curvature. This count is closely related to the number 7Cp(x) 



of prime curvatures less than x in a packing P (see Remark 2.6 1. We confirm experimentally that our heuristic 
holds for the packings Pp and Pq- We note that our heuristic does not depend on the chosen packing P - in fact, it 
yields the correct count of prime curvatures for all of the packings we checked. We summarize it in the following 
conjecture: 

Conjecture 1.2. Let Np(x) be the number of circles in a packing P of curvature less than x, and let ^fp(x) be as in 
Then as x —t 00, 

yr P (x)~L(2,z*)-Np(x) 

where L(2, £4) = 0.9159 ... is the value of the Dirichlet L-series at 2 with character Xa(p) = 1 fo r P = 1 (4) and 
X4(p) = -lforp = 3(4). 
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Samak also shows in [S2| that there are infinitely many pairs of tangent circles of prime curvature (we call these 

(2) 

kissing primes). We address the question of counting kissing primes in P via the weighted sum y/^, '(x): 

(1.4) Vp 2 \x) = £ log( fl (C))-log(a(C')) 

(c.c')es 
a{C),a{C')<x 

where S is the set of unordered pairs of tangent circles (C,C) of prime curvature in a packing P, and a(C) and 

(2) 

a(C ) denote their respective curvatures. In this case, it is less obvious what the relation is between Xffp 1 (x) and the 

number 7tp(x) of kissing prime circles in a packing P both of whose curvatures are less than x. We therefore stick 

(2) 

with \j/p ( x ) in our computation: 



(2) 

Conjecture 1.3. Let y p (x) be as in ( 1.4 ), and let Np(x) be the number of circles in a packing P of curvature less 
than x. Then 

V [ p\x)~c-L 2 (2,Xa)-N p {x), 
where Np(x) is as above and c = 1 .646 ... is given by 

These heuristics are computed by counting primes in orbits of the Apollonian group, which is possible due to 
recent results of Bourgain, Gamburd, and Sarnak in BBGSL as well as the recent asymptotic count of Kontorovich 
and Oh in |KO| for the number Np(x). Our computer experiments were conducted using Java and Matlab, and the 



programs are available at http://www.math.princeton.edu/~ksanden/ElenaKatCode.html A brief description of our 
algorithm and a discussion of its running time can be found in Section|4] 

1.1. Arithmetic structure of the Apollonian group and its orbit. Since all of the computations and claims in 
this paper concern the orbit of the Apollonian group A acting on a vector v £ Z 4 , we recall the description of the 
orbit modulo d for any integer d from [Fl |. We use this description throughout Sections|2]and[3] 

Theorem 1.4. (Fuchs): Let be an orbit of A acting on a root quadrupl^of a packing, and let 0^ be the reduction 
of this orbit modulo an integer d > 1. Let C = {v 7^ |-F(v) = 0} denote the cone without the origin, and let be 
C over 'L/d'L: 

C d = {v £ TLjdTL |v # (d), F(v) = (d)} 
Write d = d\di with (^2>6) = 1 and d\ = 2"3 m where n,m > 0. Write d\ — v\V2 where V\ — gcd(24,d\). Then 

(i) The natural projection 6& — > 0^. x is surjective. 

(ii) Let n : C c i l — > C V] be the natural projection. Then ff^ — n \& Vx ). 

(iii) The natural projection 0^ — > Ylp r \\d 2 &p r ^ s surjective and p r = C p r. 

This result is obtained by analyzing the reduction modulo d of the inverse image of the Apollonian group A in the 
spin double cover of SO p. We note that Theorem 1 1 ,4| implies that the orbit 6 of A has multiplicative structure in 



2 

A root quadruple of a packing P is essentially the 4-tuple of the curvatures of the largest four circles in P. It is well defined and its properties 
are discussed in |GLMWY|. 
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reduction modulo d = Y[p r \\d p' ar, d that it is completely characterized by its reduction mod 24, or by ^24 in our 
notation. This explains the dependence on P24 in Conjecture 

Acknowledgements: We thank Peter Sarnak, Alex Kontorovich, and Kevin Wayne for many insightful conversa- 
tions and helpful suggestions. 

2. PRIME NUMBER THEOREMS FOR ACP'S 

In JBGS I, Bourgain et.al. construct an affine linear sieve that gives lower and upper bounds for prime and 
almost-prime points in the orbits of certain groups. In this section, we use their analysis to predict precise asymp- 
totics on the number of prime curvatures less than x, as well as the number of pairs of tangent circles of prime 
curvature less than x in a given primitive Apollonian packing P. The conditions associated with the affine linear 
sieve for A are verified in [BGS]. We recall the setup below. 

Let a„ — #{circles of curvature n in a bounded packing P], and note that a n is finite since the number of circles 
of any given radius can be bounded in terms of the area of the outermost circle. We consider 1 < n < x and note 
that the sum 

Y^a„ =N P (x), 
n 

where Np(x) is the number of circles of curvature less than x and is determined by the asymptotic formula in |KO| 
(see Lemma [2~2] ). Key to obtaining our asymptotics is computing the averages of progressions mod d of curvatures 
less than x, where d > 1 ranges over positive square-free integers of suitable size. To this end, we define 

%d = a " 

n=Q(d) 

and introduce a multiplicative density function j3 (d) for which 

X d =P(d)-N P (x) + r(A,d) 

where the remainder r(A,d) is small according to the results in |BGS |. 

In the case of ACP's, we define /3 (d) as follows. Let 6 be an integral orbit of A, and let 6& be the reduction of 
& modulo d for a square-free positive integer d. Then 

(2 1) Bid)- #{yeffAV}=Q} 
{IA) PjW #{veft} 

where vj is the jth coordinate of v. We recall from Theorem |1.4| that the orbit 6& has a multiplicative structure 
which carries over to the function j3/ so that 

P \d 

Thus in order to evaluate Pj(d) for arbitrary square-free d, we need only to determine Pj(p) for p prime. This is 
summarized in the following theorem. 
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Lemma 2.1. Let d = Y[pi be the prime factorization of a square-free integer d > 1. Then 

(i) Pj(d) = UPj(Pi)forl<j<4. 

(ii) For p j^2, we have 

Pj(p)=Pk(p)forl<j,k<4. 

(iii) For any orbit & there exist two coordinates, i and j, such that 

ft(2)=/3y(2) = l, 
p k (2)=0fork^i,j. 

We say that the ith and jth coordinates are even throughout the orbit, while the other two coordinates are 
odd throughout the orbit. 

(iv) For p^2, let J3 (p) = pi(p)for 1 < i < 4. Then 

{— |-r for p = 1 mod 4 
^ forp = 3mod4 



Proof. The statements in (i) and (ii) follow from Theorem 1 .4 Let v be the root quadruple (the quadruple of the 
smallest curvatures) of the packing P. To show (iii), note that any quadruple in a primitive integral ACP consists 
of two even and two odd curvatures (see [Sn| for a discussion). Without loss of generality, assume v = (1, 1,0,0) 
mod 2, so i = 1 and j = 2 in this case. Since the Apollonian group is trivial modulo 2, we have that every vector in 
the orbit is of the form (1,1,0,0) mod 2, so we have what we want. 

To prove (iv), we use results in BF11 and recall from Theorem 1.4 that & p is the cone C p for p > 3. Thus the 
numerator of j3 (p) is 

#{v e e P I vj = o} = #{(vi , v 2 , v 3 ) eFp-{o}|F(vi, v 2 , v 3 , o) = 0} 

where F is the Descartes quadratic form and p > 3. So the numerator counts the number of nontrivial solutions to 
the ternary quadratic form obtained by setting one of the v; in the Descartes form F(\) to 0. Similarly we have that 
the denominator of j3 (p) is 

#{v € ff p ) = #{(vi , v 2 , V3, v 4 ) e ¥ p - {0} | F( Vl , V 2 , V3, v 4 ) = 0} 

where p > 3. So the denominator counts the number of nontrivial solutions to the Descartes form. The number 
of nontrivial solutions to ternary and quaternary quadratic forms over finite fields is well known (see JC)> f° r 
example). Namely, 

p 3 + p 2 — p — 1 for p = 1 mod 4 
p 3 — p 2 + p — 1 for p = 3 mod 4 

for p > 3, and 



(2.3) #{(v 1 ,v 2 ,V3,V4)eF^-{0}|F(vi,v 2 ,V3,v 4 ) = 0}= . ; , 



(2.4) #{(vi , v 2 , v 3 ) G F ; ; - {0} | F(v 1 ,v 2 ,v 3 ,0) = 0} = - 1 for all odd primes p. 



SOME EXPERIMENTS WITH INTEGRAL APOLLONIAN CIRCLE PACKINGS 



7 



Combining (2.3 i and (2.4 1, we obtain the expression in (2.2 1 for p > 3. For p = 3, we compute G p explicitly and 
find that there are two possible orbits of A modulo 3 which are illustrated via finite graphs in Fig. [3] and Fig. |4] 
Both of these orbits consist of 10 vectors v G Z 4 . In both orbits, 4 of the vectors v have v; = for any 1 < i < 4. 
Thus J3(3) = I as desired. □ 




Figure 3. Orbit I modulo 3 




Figure 4. Orbit II modulo 3 



In the following two sections, we use this setup to produce a precise heuristic for the number of circles of prime 
curvature as well as the number of pairs of tangent circles of prime curvature less than x in a given ACP. 



2.1. Predicting the prime number theorem for ACP's. In order to compute the number of prime curvatures in 
an ACP as proposed in Conjecture 1 .2 we use the setup above paired with properties of the Moebius function to 
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pick out primes in the orbit of A (see (2.7 1). We use the asymptotic in [K0| for the number Np(x) of curvatures 
less than x in a given packing P: 

Theorem 2.2. (Kontorovich, Oh): Given a bounded Apollonian circle packing P, there exists a constant cp > 
which depends on the packing, such that as x — > oo, 

Np(x) ~cp-x 5 , 

where 8 — 1.30568 . . . is the Haussdorf dimension of the limit set of A acting on hyperbolic space. 

For the purpose of our computations, we will need a slightly stronger statement of Theorem |2.2| as we will 
sum over each coordinate of the points in the orbit of A separately. Namely, each circle in the packing is uniquely 
represented in the orbit as a maximal coordinate of a vector v in Z 4 . We would like to know how many circles 
there are of curvature less than x that are represented in this way in the z'th coordinate of a vector in the orbit. We 
denote this by N P (x) : 

(2.5) 



veff 
v*<x 

where v* denotes the rth coordinate of v G Z 4 which is also a maximal coordinate of \[^] To this end we have 
Lemma 2.3. Let N P '\x) and Np(x) be as above. Then 
(2.6) 

as x approaches infinity. 



4%) ~ N® (x) ~ N® (x) ~ N P 4) (x) ~ 



Proof. The computation in IIKOI of the main term in the asymptotics in Theorem 2.2 relies on the Patterson- 
Sullivan measure on the limit set of the Apollonian group A. In order to prove Lemma |2~3] we show that this 
measure is invariant under transformations on the coordinates of a vector v in an orbit G of A. 

To this end, let G be the group of permutations of the coordinates vj,...,V4 of a vector v € 6. The group G is 
finite and its elements can be realized as 4 by 4 integer matrices. For example, the matrix 

/ 1 \ 

10 

10 

\ 1 j 

switches the first and second coordinates by left action on v T . Let 

L=(G,A) 

be the group of 4 by 4 matrices generated by the Apollonian group A together with G, and note that each element of 
G normalizes A. For example, if M is as above, we have M~ l S\M — S2, where Si and S2 are as in ( 1.2 1. Similarly, 



M 



It is possible that there is more than one i for which the i'th coordinate is maximal 
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any element of G switching the zth and /'th coordinate of v conjugates Sj to Sj in this way. Thus L/A is finite, 
and so L is a finite extension of A. In particular, this implies that the Patterson-Sullivan measure for L is the same 
as for A. Since L is precisely an extension of A by the permutations of the coordinates of v, we have that the 
Patterson-Sullivan measure is invariant under G. Together with Theorem 2.2 and its proof in |KO], this proves the 
lemma. □ 



Since we are interested in counting the points in 6 for which v* is prime, we sum over points for which v* is 
modulo some square-free d. It is convenient to count primes in the orbit of A with a logarithmic weight. To this 
end, we consider the function 

f log/? if n = p l 
A(n) = < 

I otherwise 

for which it is well known that 



(2.7) 



A(») 



d\n 



where fi (d) is the Moebius function. Using this we write down a concrete expression for the number of prime 
curvatures less than x in a packing P counted with a logarithmic weight: 



Lemma 2.4. Let v* be the ith coordinate of a vector v in & such that v* is the maximal coordinate of v, and let 
y/>(x) be as before. Then 



(2.8) 



£ £A(v?)+0(*) 



l<i<4 \e0 



The sum in (2.8 1 is a count of all circles whose curvatures are powers of primes. Including powers of primes in 
our count will not affect the final answer significantly. Namely, let Np- 1 (x) be the number of circles in a packing P 
whose curvatures are less than x and perfect squares. Note that 

N^(x) =0(jc). 



This is insignificant compared to the count of all curvatures in Theorem 2.2 so the sum in (2.8 i is the correct one 
to consider. Denote by D < Np(x) be the level distribution from the analysis in [BGS | - i.e. our moduli d are taken 



to be less than D. We combine Lemma |Z4| Lemma 2.3 and the expression for A(n) in (2.7 1 to get 



Wp(x) 



£ £ 2>(«0iog</+o(x) 



1<K4 ve^ d\v\ 
v*<x 



(2.9) 



£ £ 5>Miog«/Ei- E E J>(<01ob«*£i+o(*) 

1<;<4 vtfJ d<D v*=Q(d) KK4 ve^ d>D vf=0 (d) 

v* <x 



Assuming that the Moebius function fl(d) above becomes random as d grows, the sum over d > D in (2.9 1 is 



negligible, and we omit it below. We proceed by rewriting the sum over d < D in (2.9 1 using the density function 
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[5(d) in (pTl. Recall that the analysis in MBGS1 and BKOl gives us that 



£ a n = P{d) -c P x S +r(A,d) 

n=0(d) 

where r(A,d) is small on average. In particular, 

Y j r(A,d)=0(x 8 -^) 

d<D 

for some £o > 0. Paired with the assumption that jj. is random, this evaluation of the remainder term allows us to 
rewrite (12. 9b as follows: 



N' P (x) 



£ ^{dMd)\ogd+o(x 5 -^) 



l—i A 

1<K4 ^ d<D 

(2.10) = - Np ^ £ £ ft^Wlogd + Otjt*- 1 *) 

4 1</<4c/<D 

To compute the innermost sum in the final expression above, note that 

(2.11) £ pi{d)n(d)logd = £ fr(d)ii{d)\ogd- £ A'(^)MWlogrf. 

t/<£> rf>0 

Assuming once again that the sum over d > D is insignificant due to the conjectured randomness of the Moebius 
function, we have that the sum over d < D in ( |2.1 1[ ) can be approximated by the sum over all d. With this in mind, 



the following lemma yields the heuristic in Conjecture 1.2 



Lemma 2.5. Let j5j(d) be as before. We have 



£ £pi(d)n(d)logd = 4-L(2, X4 ) 

l<i<4d>Q 



where L(2,Xa) = 0.91597 . . . is the value of the Dirichlet L-function at 2 with character 

Xa{p) = 



1 if p = 1 (mod 4) 
-1 if p = 3 (mod 4) 



Proof. We introduce a function 

f(s)=Y,Pi(dMd)d- s , 

d 

and note that its derivative at is precisely what we want: 

f(0) = -Y t Pi(d)n(d)logd. 

d 
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Since the functions j3, ji, and d s are all multiplicative, we may rewrite f{s) as an Euler product and obtain 

f{s) = U(l-P,(p)p-') 

p 

= nc-p-') i^f 

P I 

p 1 r 

= Z-\s+\).H{s) t 



where H{s) = Y\p{^ ~ fii{p)p s )(l ~ P s ) 1 is holomorphic in 9t(i) > 1/2. Differentiating, we obtain 

/'(o) = -c'(i)r 2 (i)-^(o)+r 1 (i)-// , (o)=//(o) 

since — £'(1)£~ 2 (1) = 1 and £ _1 (1) =0. Thus it remains to compute 



^)-riT— r 



Part (iii) of Lemma 2.1 says j8,(2) = j8/(2) = 1 for two coordinates 1 < i,j < 4. For these two coordinates 1 



Pi(2) = and so H(0) = 0. Otherwise j3,-(2) = and we have 

2 2 

11 p 2_l 11 „2 + 1 

= 2-L{2,Xa)- 

Thus the sum we wish to compute is 4 • L(2, £4), as desired. □ 



Lemma [23] implies that the contribution of the two of the coordinates that are even throughout the orbit to the 
sum in ( |2.10| i is 0, and the contribution for the other two coordinates is 

^■4-L(2,X4)=N P (x)-L(2,X4), 



yielding the predicted result in Conjecture 1.2 



Remark 2.6. It is well known that Kp{x) ~ as x — > °°. Thus Conjecture 
7t P (x): 

L(2,X4)-N P (x) 



1.2 



can also be stated in terms of 



Tlp(x) 



\ogx 
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Figure 5 . Prime Number Heuristic for P c 
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Figure 6. Prime Number Heuristic for P B 



Since our computations rely on the multiplicativity inherent to the reduction mod d of the Apollonian group A, 
our heuristic is independent of the chosen packing P in which we count prime curvatures. This is confirmed by our 
data: Fig. [5] and Fig. [6] above show the graphs of 



(2.12) 



Np(x) 
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where P = Pq and Pg as in Section [T] and x < 10 8 . In both cases, |Z12J) converges to y — L(2,%a) as predicted in 



Conjecture 1.2 



2.2. Predicting a prime number theorem for kissing primes. In this section, we use the analysis in MBGSL as 



well as the conjectured randomness of the Moebius function to arrive at the heuristic in Conjecture 1.3 for the 
number of kissing primes, i.e. pairs of tangent circles both of prime curvature less than x. With the same notation 



as in Section 2.1 we would now like to count the points in & for which v* and v; are prime for some j ^ i, so we 
sum over points for which either v* or V/ is modulo some square-free d. To do this we need the total number 
of pairs of mutually tangent circles of curvature less than x in a packing P. If Np(x) is the number of circles of 
curvature up to x as specified by Theorem |2.2| it is not difficult to see that 

(2.13) #{pairs of mutually tangent circles of curvature a < x} = 3 • Np(x) 

since each circle of curvature a in the packing is tangent to a distinct triple of circles of curvature < a. This can 
be seen via a simple proof by induction. We again employ the function A(n) in order to write down a concrete 
expression for the number of kissing primes less than x in a packing P. 



Lemma 2.7. Let v* and Vj be two distinct coordinates of a vector v in 6", where v* denotes the maximum coordinate 

(2) 

of\, and let \ff p (x) be as before. Then 



(2.14) 



V?\*) = EJ>(v?)A(vy) + 0(x). 



V£ff jyti 



Again, the sum in (2.14i is a count of all mutually tangent pairs of circles whose curvatures are powers of 
primes, but by a similar argument to that in Section[2~T|we have that including powers of primes in our count does 



not affect the final answer significantly. Note that in order to evaluate ( 2. 14 1 we introduce in ( 2. 1 8 1 a function which 



counts points in for which two of the coordinates are modulo p. Denote by D < x be the level distribution 
from the analysis in [BGS | - i.e. the moduli d > 1 in the computations below may be taken to be less than D. We 
rewrite the expression in ( |2.14| ) using ( |2.7) > and get 



(2.15) 
(2.16) 
where 



E E E M(*)lDgdi E li(.dj)]ogdj 

l<;,j<4 veff \di\v* dj\vj j 

¥j v l< x 

E E(^+ E+ )+°w 

l<ij'<4 veff 
¥j V?<X 



dj<D v?=0(di) dj<D Vj=0{dj) 
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and 



(*i>D v*=0(c/,) f/j>D Vj=0{dj) 



As in Section 2.1 we omit E + in (2.15 i under the assumption that ji behaves randomly for large values of d. Along 



with the results about the remainder term in the sieve in [BGS|, the expression in (2.15 1 becomes 



(2.17) Np(x) £ £ ft( ? ^)p/ 7 ^) g ((4,^))/i(di)M(^)log£/aog^+0(^-*), 

l<i,/<4[rfi,d;]<D' \\ a ii a })J \\ a U a j)J 
¥j 

where f5j(d) is as before, [di,dj\ is the least common multiple of di and dj, and (di,dj) is their gcd. The function g 
above is the ratio 

#{v G {dudj) |vj = 0((<My)) and =0((4rf / ))} 



(2.18) 



g{{di,dj)) 



#{ve^)} 

where (di,dj) is square-free in our case. Note that g(c/) is multiplicative outside of the primes 2 and 3 by Theo- 
rem |1.4| so 

g(d) = UP 

p\d 

and we need only to compute g(p) for p prime in evaluating the sum above. 



Lemma 2.8. Let g(j>) be as before where p is a prime. Then 
1 if both vi and Vj are even 
if at least one of Vi or vj is odd 



(0 8(2) - 
(ii) 8(P) 



(p+i) 2 



p 2 +i 



for p = 1 mod 4 
for p = 3 mod 4 



Proof. To prove (ii), we note that Theorem jl ,4| implies that the numerator of g(p) is 

#{v eff d \v!=0 (d) and v 2 = 0{d)}=#{ (vj ,v 2 )e¥ z p - {0} | F(vi , v 2 , 0, 0) = 0} 

for p > 3. Thus it is the number of non-trivial solutions to a binary quadratic form with determinant (the Descartes 
form in ( |1.1[ ) with two of the v/, v 7 set to 0), and so 

#{v € ^p|vi = mod p and V2 = mod p} = p — 1 

for all p > 3 (see (C |, for example). In the case p — 3, we observe that in both of the possible orbits of A mod 3 in 
Figures |5] and ^ we have g(3) — jq as desired. Part (i) follows from the structure of the orbit G 2 as observed in 
the proof of Lemma 2.1 □ 



Denote by v(dj,dj) the gcd of di and dj (we write just v from now on and keep in mind that v depends on d{ 
and dj). Note that 



[di,dj] 
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where and are relatively prime, so we rewrite the expression in (2.17 i below: 

N P (x 



(2.19) 



, E E p i ( d Ap ] ( d ±)g{v) i i{d i ) i i{d j )\ogd i \ogd ] . 

* \<i.j<4[dj.dj]<D' V v / V v 



To compute the sum in (2.19 1, we use a similar argument to that in Section 2.1 We note that the inner sum in the 



expression above is equal to 
(2.20) 

E E Pi Pj givMdiMd^iogdiiogdj - £ Pt (%) pj (^) givMdiMd^iogdiiogdj 

di>odj>o v v / v v / \di,dj\>iy \ J \ J 

where we assume the sum over [di,dj] > D' is insignificant by the conjectured randomness of the Moebius function, 



and thus the sum over all di and dj is a good heuristic for the sum in (2.19 1. We compute the infinite sum in the 



following lemma and obtain the heuristic in Conjecture 1.3 



Lemma 2.9. Let Pi(d), Pj(d), and g(v) be as before. We have 

E E E Pi (if) Pi (if) g(v)n(di)n{dj)\ogdilogdj = 8-L 2 (2,£ 4 ) ■ U (l 

l<i,j<*di>0dj>0 V V / V V / P=3(4) 



Proof. We introduce the function 



f(si,sj) = (f) Pj ( d i ) givMdiMdjWtf, 

dj dj \ ' \ 



and note that 



d 2 f(Si,Sj) 



dsjdsj 



= LLP (^)p( d i)g{v)li{di)li{d ] )\ogd i \ogd h 
(o.o) ^ dj V v / V v 



which is a good heuristic for the sum in (2.19i as x tends to infinity. The difficulty in computing this is the 



interaction of dj and dj in g(v). To this end, write dj = ve; and dj = Ve ; - where = 1. This gives us the 

following formula for f(si,Sj). 

f( Si , Sj ) = J>(v) E PiiedPjiejMyeiMvejXve^iyejYJ 
(2.21) = ^g(v)v-M£ £ rtmMetfjiejMveiMvejWe 1 / 

v ei,ej m \{e b ej) 

= £g(v)v-( Ji+ ^£/i(«)M 2 (vm)m-( Si+s ^ 



£ PiimbJPjimbjMbiMbjWb*; 

(&y,Vm)=l 
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where 



A(v,m, Si , Sj ) = £ PiibOpjibjMbiMbjWb/' 

{bj,Vm) = \ 

(bj,vm) = l 

= n ((i-AW)(i-ft-0'>7 ,7 ))~ 1 

p ; |vm 

(2-22) ■ ri( 1 -AW)(l-ft(Py)/'7 , 

. r . (Sl+1)n lz^ r , (S;+1)n lz£d£M! 

Pi ' P; 1 



where 



We write 



which gives us 



(2.23) 



Py(si,sj) = J! ((1 - Pi(Pi)PT")(l -PjMpJ"))' 1 ■ 



Pi\ v 
Pj\ v 



P m (si,sj) = n ((i - A -Pj(pj)p7 j )) 



-1 



Pil' 
Pj\" 



B (*.^) = I1— — =S=i~Il 



C(v, Si , Sj )= 8 (v)v-^ +s J\ 
D(m,Si,Sj) = v(m)j3i(rn)fij(m)m~( Si+s ^ 



v (m,v)=l 



Note that the sums over v and m, as well as B(si,Sj) converge and are holomorphic. We now compute the desired 
derivative. Write 

Gij(y,m,Si,Sj)=Y,C(v,Si,Sj)-P v (si,Sj) £ D(m,Si,Sj)-P m (si,Sj). 

v (m,v)=l 
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Then we have 

d 2 f( Si ,Sj) 



(2.24) 



r 2 (i) 



dsidsj 



(0.0) 



d 2 (Gij(v,m,Si,Sj)) 



Gi j(v,m, 0,0) -.8(0,0) 



(0,0) 



C'(l) 

C 2 (i) 



Since C~ 2 (l) = 0, we need not compute the partial derivative of G. If j8,(2) = ]3/(2) = in the expression for B, 
we have that B(0,0) = 4-L 2 (2,^ 4 ), since = 2 -L(2,^ 4 ) as in Section j2~lj This holds for the sums in (2.19 1 
over (v*, Vj) and (v*, v;) where the jth and jth coordinates in our orbit are everywhere odd. 



The contribution to (2.19 1 from the terms where v* or v, is even is 0: 



Lemma 2.10. Let & be an orbit of the Apollonian group. Given that the i-th or j-th coordinate of each vector 
v £ & is even, we have 

G Lj {v, 111,0,0) -B(0,0) =0 



Proof. Recall that two of the coordinates of the vectors in & are always even, and two are odd. We write 
Gij(v,m, 0,0) -5(0,0) - £(0,0)(£C(v,0,0)-P v (0,0) £ D(m,0,0) -P m (0,0) 

\2|v (m,v)=l 

(2.25) + £ C(V,0,0)-P V (0,0) £ D(/n,0,0)-P m (0,0) ) . 

(v,2) = l (m,v) = l / 



Case 1: only one of the i-th and j-th coordinates is odd throughout the orbit. 



Recall from Lemma 2.8 that g(2) =0 in the case that only one of the coordinates (z',j) is even throughout the 
orbit, and g(2) = 1 if both coordinates are even throughout the orbit. So, if only one of the coordinates (i,j) is 
even throughout the orbit, we have 

B(0,0)-[£C(v,0,0)-Pv(0,0) £ D(m,0,0)-P m (0,0) J =0. 

\2|v (m,v) = l ) 
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Recalling that j3,-(/?) = fij(p) for p > 2 from Lemma 2. 1 we have 
(2.26) 



5(0,0) £ C(v,0,0)-P v (0,0)( £ D(m,0,0)-P m (0,0) 



v,(v,2)=l 
( (1-Pt(pi))(l-Pj(pj)) 



(m,v) = l 



^ C(v,0,0)-Pv(0,0)( X O(m,0,0)n(l-ft(/>)r 2 

v(v,2)=l 



(m,v)=l 
(m,2)=l 



+ 5(0,0)- £ C(v,0,0)-P v (0,0)( ^ p(2m)ft(2m)Pj(2m)-P m (p,0) 

\(v,2) = l (2m,v)=l 

= + 
= 



since Lemma 2. 1 implies either 1 — /3j(2) = or 1 — j3;(2) = in the first term, and either /3, (2m) = or J3y (2m) = 



in the second term in (2.26 1. We now compute the expression in (2.25 i in the case that both the ith and jth 
coordinate are even throughout the orbit. 



Case 2: the i-th and j-th coordinates are both even throughout the orbit. 



In this case, Lemma 2.8 implies that g(2v) = g(v) and that C(2v,0,0) = C(v,0,0) for odd v. Also, we again 



have that pi(p) — Pj{p) for p > 2, so the first sum in (2.25 1 is 
(2.27) 



5(0,0)- [£C(v,0,0)-P v (0,0)( £ D(m,0,0)-P m (0,0) 



\2|v 



(m,v)=l 



£ (C(2v,0,0)n(l-ftb)r 2 )( I (D(m,0,0)n(l-ft(p))- 

\(v.2) = l p|v (m,2v) = l p|m / 

( £ (C(v,0,0)n(l-A(rt)- 2 )( E (D(m,0,0)n(l-A(^))- 2 )))-i 2 (2,Z4) 

\(v,2) = l p|v (m,2v) = l p|m / 



To compute the second sum in (2.25 1, we note that /5,-(2m) = j3y(2m) = /3,(m) by Lemma 2.1 and write 



5(0,0) £ C(v,0,0)-P v (0,0)( £ D(m,0,0)-P m (0,0) 

\(v,2)=l 

M 1 (v,0,0)+M 2 (v,0,0), 



(m,v)=l 
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where 



M 1 (v,0,0)= £ C(v,0,0)-Pv(0,0)( £ }i{m)pfim)-Yl{\-fr{p))- 2 ) -5(0,0) 

\(v,2) = l (m,2v) = l p|m / 

M 2 (v,0,0)=( £ C(v,0,0)-P v (0,0)( £ M(2m)A 2 (m)-n(l-A(p))" 2 ) J ^'(0,0) 

\(v,2) = l (m,2v)=l pjm / 



where 

Note that, since 1 -$(2) = 1 ~Pj(2) = 0, we have fi(0,0) = in this case, so 

Mi (v, 0,0) =0. 

On the other hand, 
(2.28) M 2 (v,0,0) 

= ( £ C(v,0,0)-P v (0,0)( £ -M(m)A 2 (m)-n(l-AWr 2 )]-B'(0,0) 

\(v,2) = l (m,2v) = l p\m / 

= -( I (C(v,0,0)n(l-A(p))- 2 )( £ (Z)(m,0,0)n(l-A(p))- 2 ))yi 2 (2,Z4) 

\(v,2)=l p|v (m,2v) = l p\m ) 

Combining ( |2.27[ i and ( |2.28| l, we have that the expression in ( |2.25| > is 0, as desired. □ 



Therefore, we have that in the contributing terms of (2.19 1 both v* and v ; are odd. This makes up two of the 



terms in ( |2.19| l, so we have that the sum we wish to compute is equal to 

(2.29) 2 • Gy(v,m, 0, 0) • 4 • L 2 (2, X a) = 8 • L 2 {2, X a) ■ Gij{v,m, 0, 0) . 

It remains to compute Gjj(v,m,0,0) in this case. Recall that j3/(2) = jB/(2) = in this case, and that j3,(p) = j3 y - (/? ) 
for p > 2. We have 

£ jU(m)/3 2 (m)]~[(l-/3,-(/?)r 2 

(m,v)=l p|m 

- m- Airt- ft2(p) 



(i-AW^V (i-A(p)) 2 



n i-- n i-^-nfi- m ^ ' 



We write 



„ n 1 1 n 1 (P + 1 ) 2 

° = 11 !-72 11 l ~ (p 2_ p) 2 
p=l(4) ^ p=3(4) ^ ^ 



20 

and get 



ELENA FUCHS AND KATHERINE SANDEN 



(2.31) 



G (J (v,m,0,0) 

,.£c(v,o,o,, v( o,o,.n(.-^ ? ) 



-1 



g(p) 



a -in (1 _ j8 2 (/ , )(1 _ j8(/ , )) -2 )(1 _ j8(/ , ))2 

rr, . gip) 

V + (l-j32(p)(l- j 3(p))-2)(l-^(p))2 



1 



-1 



p 2 + l 



p=l(4) V ^ / p=3(4) 



Therefore our infinite sum is equal to 
(2.32) 



p 2 + l 



//-2p 3 -2p-l 



p=3(4) A J 



as desired. 



□ 



Conjecture 1.3 follows from our assumption that the Moebius function is random and from Lemma ??. namely, 
we predict 

-c-L 2 (2, X a), 



v4 2) w-8.^.^(2, Z4) . n i 2 



p=3(4) 



P {p -If 



where c = 1 .646 ... as in Conjecture 1 .3 



_V p 2) (x) / N< 2) (x) 



n0.49 

-0.48 
-0.47 
0.46 
0.45 
0.44 
0.43 
42 



1C 

x10 6 



FIGURE 7. Prime Number Theorem for Kissing Primes for the packing Pg 
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_i|i<?(x)/N<?>(x) 



0.48 



-0.47 



0.46 



4 6 
x 



-0.45 



-0.44 



'.43 



x 10 



Figure 8 . Prime Number Theorem for Kissing Primes for the packing Pq 



As with our heuristic for the number of prime curvatures less than x in a packing, this count does not depend on 
the packing P. Fig.[7]and Fig.[8]above show graphs of 

for P = Pb and Pc, and x < 10 8 . Both of these graphs converge to a = ' — l 2 '* 4 ^ as conjectured, even though they 



converge slower than in the case of counting circles of prime curvature in Section 2.1 



3. Local to Global Principle for ACP's 

In this section we present numerical evidence in support of Conjecture |1.1| which predicts a local to global 
principle for the curvatures in a given integral ACP Since the Apollonian group A is small - it is of infinite index 
in Of(Z) - it is perhaps surprising that its orbit should eventually cover all of the integers outside of the local 



obstruction mod 24 as specified in Theorem 1 .4 Proving this rigorously, however, appears to be very difficult. An 
analogous problem over Z would be to show that all large integers satisfying certain local conditions are represented 
by a general ternary quadratic form - this analogy is realized by fixing one of the curvatures in Descartes' form and 
solving the problem for the resulting ternary form. While this problem has been recently resolved in general in MCol 
and IIDSI , even there the local to global principle comes in a much more complicated form, relying on congruence 
obstructions specified in the spin double cover. Our conjecture, which therefore has the flavor of Hilbert's 11th 
problem for an indefinite form, predicts a local to global principle of a more straightforward nature. 

Our computations suggest that this conjecture is true and we predict the value Xp in the examples we check. 
We consider the packings Pg and Pc introduced in Section[T] Recall that Pg corresponds to the orbit of A acting on 
(—1,2,2,3), and Pq corresponds to the orbit of A acting on (—1 1,21,24,28). 
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In order to explain the data we obtain in both cases, we use Theorem |1.4| to determine the congruence classes 



(mod 24) in the given packing. Recall that the Apollonian group A is generated by the four generators S,- in ( 1.2 1. 
We can view an orbit of A modulo 24 as a finite graph ^24 in which each vertex corresponds to a distinct (mod 24) 
quadruple of curvatures, and two vertices v and V are joined by an edge iff S,-v = V for some 1 < i < 4. Recall 



from Theorem 1.4 that for any orbit & of the Apollonian group, 



(3.1) ^24 = ^X^3, 

so the graph ^24 is completely determined by the structure of 1^3 and 0%. There are only two possible orbits 
modulo 3, pictured in Fig.|3]and Fig. |4] There are many more possible orbits modulo 8, and we provide the graphs 
for these orbits in the case of Pg and Pc in Fig. [9] 




Figure 9. Orbits of P B and P c modulo 8 



Note that each vertex in Sfs and ^3 is connected to its neighboring vertices via all of the generators S 1 ,-. Therefore 
the curvatures of circles in a packing modulo 24 are equally distributed among the coordinates of the vertices 



in ^24- Combined with Theorem 1.4 this lets us compute the ratio of curvatures in a packing which fall into a 
specific congruence class modulo 24. Namely, let ^24 (P) be the orbit mod 24 corresponding to a given packing 
P. For w £ &2A (P) let Wj be the ith coordinate of w. We define y(n,P) as the proportion of coordinates in &ii\(P) 
congruent to n modulo 24. That is, 

I t U # {w£&2A(P)\w i = n} 



(3.2) r(n,P) = 



4 ■#{ we ^24 OP)} 
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With this notation, a packing P contains a circle of curvature congruent to n modulo 24 iff y(n,P) > 0. Given (3.1 
we express y as follows: 

(3 3) yinP) = Et.*{"e*24(P)k,- = n} 

lti # {w E t \ Wj = n (3)} -#{w E0 3 \w,=n (8)} 

4-#{we &%\ -#{we 

The significance of 7 in the case of any packing (not only the two we consider) is explained in the following lemma. 

Lemma 3.1. Let Np(x) be as before, let C be a circle in an integral Appolonian packing P and let a(C) be the 
curvature ofC. Then 

£ l~y(n,P)'Np(x) 

CeP 

a(C)<x 

a(C)=n (24) 



This follows from Theorem 1.4 Note that in general the orbits @%{P) and ffj,(P) have 4 and, respectively, 10 
vertices in the corresponding finite graph^] Therefore ^24 always has 40 vertices, and the ratio in (3.3 1 is easily 
computed using this graph. With this in mind, we observe the following about the packing Pg- 

Lemma 3.2. Let Pg,24 denote the possible congruence classes of curvatures mod 24 in the packing Pg, and let 
Np B (x) be as in Theorem . 



2.2 



Then we have 



(i) N Pb (x) - c Pb -x 8 , where c Pb = 0.402 . 

(ii) Pbm = {2,3,6, 11, 14, 15, 18,23} 



(iii) 



y(2,P B ) = 4 y(l4,P B ) ;1 



y(3,P B ) = ^ rti5,ft) = i 

7(6,^) = ^ y(18,/fc) = i 

7(11,^) = 4 7(23,P B ) = 4 



(iv) For 10 6 < x < 5 ■ 10 s , let X24 denote x mod 24. If X24 G -Pb.24 then x is a curvature in the packing Pg. 
Part (iv) is an observation based solely on our computations using the algorithm described in Section [4]- these 



are illustrated in the histograms below. The first three parts follow from computations combined with Theorem 2.2 



and Lemma 3.1 Note that y(n, Pb) = y(n + 12, Pb) - for this particular packing, one can hence express the local 
obstructions modulo 12 rather than modulo 24. Whenever this is the case for an integral ACP, we will find that there 

4 There are only two possible orbits mod 3, but many more mod 8. We examine just two such orbits here. 
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are eight congruence classes modulo 24 in the curvatures of the circles. Graham et.al. observe this in IGL MWYl 
where they compute which integers less than 10 6 are "exceptions" for Pg~ they find integers which satisfy these 
local conditions for Pb modulo 12 but do not occur as curvatures in the packing Our data extends their findings 
(we consider integers up to 5 x 10 8 ), and shows that all integers between 10 6 and 10 7 belonging to one of the 
congruence classes in part (ii) of Observation [X2] appear as curvatures in the packing Pb- 

The following histograms illustrate the distribution of the frequencies with which each integer in the given 
range satisfying the specified congruence condition occurs as a curvature in the packing Pb- The frequencies seem 



normally distributed, and the means of the distributions can be computed, as we do in ( 3.6 1. The variance, however, 
is much more difficult to predict at this time - explaining the behavior of the variance as we consider larger integers 
would shed more light on our local to global conjecture. Note that there are no exceptions to the local to global 
principle in this range whenever is not a frequency represented in the histogram (i.e. each integer occurs at least 
once). There are several other frequencies not represented in the histograms for both Pb and Pq (these show up as 
gaps in the graphs), and an explanation of this aspect would be interesting. 




1000 



m . 2 (mod 24); 10 6 s m < 10 8 m . 3 (mod 24), 10 6 < m < 10 1 



There is a small error in the computations of Graham et.al. - we have found that the integer 13806 = 6 (12) does not appear as a curvature in 
Pb- Their results do not reflect this. 
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Variance: 1 2536 
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lean: 407.32 
ariance: 12536 
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Mean: 271.64 
Variance: 6083.4 
Number of exceptions: 
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Mean: 271.55 
Variance: 6133.9 
Number of exceptions: 



400 600 
m. 18 (mod 24), 10 s sm< 10 8 
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Mean: 406.92 
Variance: 12604 
Number of exceptions: 



700 800 900 



m . 23 (mod 24), 10° s m < 10° 

We do the same analysis for the packing Pq. In this case, we must consider much larger integers than in the 
case of Pb in order to get comparable results. This can be explained partially by the fact that the constant cp in 
Kontorovich-Oh's formula 



(3.4) 



Np{x)~c p -x s 
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is much smaller for the packing Pq than the packing Pp since the initial four circles in Pq are much larger than 



the initial four in Pb (See part (i) of Lemmata 3.2 and 3.3 I. Specifically, cp c = 0.0176. . . and cp B — 0.402 . 
However, our data suggests that the proposed local to global principle should hold for this packing as well. 



Lemma 3.3. Let Pq.ia denote the possible congruence classes mod 24 in the packing Pq, and let Np c (x) be as in 
Theorem\2.2\ Then we have 



(i) Np c (x) ~ cp c ■ x s , where cp c = 0.0176 . . . 

(ii) ^.24 = {0,4, 12, 13, 16,21} 



(iii) 



Y(0,Pc) = ^ 7(13,^)= 3 



10 " ' w 10 

7(4^) = ^ Y(16,Pc) = ^ 

Y(n,P c ) = ± 7(21,Pc) = 4 



(iv) For 10 8 < x < 5 ■ 10 8 , let X24 denote x mod 24. If X24 — 13 or X24 = 21, then x is a curvature in the packing 
Pc- 



Again, note that part (iv) is an observation based solely on our computations, while the first three parts in 
Lemma |3~3"1 rely on Lemma [3~T| and Theorem |2.2| The histograms below illustrate the distribution of the frequencies 
with which each integer in the given range satisfying the specified congruence condition occurs as a curvature in 
the packing Pc- Note that, as with Pp, the frequencies with which integers are represented in the packing seem to 
have a normal distribution. However, since the mean of this distribution is much smaller for Pc than for Pp, we 
find that is often a frequency represented in the histograms, and so there are still some exceptions to the proposed 
local to global principle in the range we consider. 




Mean: 24.31 
Variance: 40.98 
Number of exceptions: 538 



Mean: 36.47 

Variance: 69.73 

Number of exceptions: 394 
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m.0(mod 24), 4*1 s 
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m <5*10 8 
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40 60 
m. 4 (mod 24), 4*1 B 



80 

m <5*10 8 
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Mean: 24.32 

Variance: 24.35 

Number of exceptions: 536 



Mean: 72.95 
Variance: 196.97 
Number of exceptions: 



10 20 30 40 50 60 70 80 90 
m . 12 (mod 24), 4*10° s m<5*10 8 




2r 
18- 
1.6- 
1.4 
1.2- 

1 - 
0.8- 
0.6 
0.4 
0.2 



Mean: 36.49 

Variance: 71 .14 

Number of exceptions: 394 



40 60 80 100 120 140 160 
m . 1 3 (mod 24), 4*1 s < m < 5*1 8 



Mean: 48.64 
Variance: 105.32 
Number of exceptions: 



20 40 60 80 

m . 16 (mod 24), 4*10 8 < m <5*10 8 




20 40 60 80 100 

m.21 (mod 24), 4*10 8 s m <5*10 8 



As we mentioned before, the mean in each of these histograms is easily computable: let C denote a circle in an 
Apollonian packing P and let a (C) denote the curvature of C. Let / = [k, k + K] be an interval of length K and let 
x £ I be an integer. Let 

v(x)=#{C eP\a(C)=x} 
be the number of times x is a curvature of a circle in P. For and integer m>0, let 

8{m,n) =#{x € I\x = n(2A), v(x) =m}. 



^ v(x) = ^ 8(m,ri) -m. 

xeJ m>0 
x=n (24) 



Then by Lemma [3TT 
(3.5) 



The equivalence of the two sums above is easy to observe - one counts the same set of curvatures, but partitions 



them differently. In particular, the expression in (3.5 1 allows us to determine the mean of the distributions in the 



histograms above. Namely, denote by x G / an integer in some interval / = [k,k + K] of length K. Let 1 < n < 24, 
and let ji(n,P) denote the mean of the number of times x = n mod 24 is represented as a curvature in the packing 



P. Note that there are precisely K/24 integers congruent to n mod 24 in the interval /. Combined with (3.5 1, this 
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gives US 

(3.6) ^ p) ^. 7{ n,P). {M k + K)- Np{ k))_ 

K 

This formula predicts the following values for the means in P B in the range [ 1 6 , 1 8 ) ; and Pc in the range [4 x 
10 8 ,5-10 8 ): 

H(2,Pb) =h(U,Pb) =h(U,Pb) =h(U,Pb) = h(23,P b ) = 406.70. . . 
H(3,P B ) = n(6,P B ) = n(l5,P B ) = M(18,ft) =271.13... 
H(0,P C )=H(12,P C ) = 2435... 
H{4,P C ) = n{l6,P c ) =36.52... 
M(13,P C ) =73.05... 
H(2l,P c ) =48.70... 

which coincides with the means observed in the histograms. This clarifies why the mean is small for packings 
where the constant cp in the formula Np(x) ~ cp-x s is small, and why one needs to consider very large integers to 
see that the local to global principle for such packings should hold. 

This analysis can be carried out for any ACP, and will likely yield similar results. In the direction of proving 
Conjee ture |l.l| one might investigate how Xp depends on the given packing - can it perhaps be expressed in terms 
of the constant cp in ( 3.4 1? One might also ask how the variance of the distributions above depend on the packing, 
and how it changes with the size of the integers we consider - answering this would give further insight into the 
local to global correspondence for curvatures in integer ACP's. 



4. A DESCRIPTION OF OUR ALGORITHM AND ITS RUNNING TIME 



We represent an ACP by a tree of quadruples. Fig. 10 shows the first two generations of the tree corresponding 
to Pc. To generate all curvatures of magnitude less than x, we use a LIFO (last-in-first-out) stack to generate and 
prune this tree. The algorithm is as follows: 

1 ) Push the root quadruple onto the stack. 

2) Until the stack is empty, perform an iterative process: 

a ) Pop a quadruple off of the stack and generate its children. 

b) For each child, if the new curvature created (i.e. the maximum entry of the quadruple) is less than x, 
then push the child onto the stack. 



By pushing a quadruple onto the stack only if its maximum entry is less than x, we effectively prune the the tree. 
Since we know each quadruple has a larger maximum entry than its parent, we use step 2b) to avoid generating 
branches whose quadruples are known to have entries greater than x. 

Although we use the concept of a tree to generate curvatures, we note that the entire tree structure is not 
necessary to store such curvatures. Instead, we store the curvatures in a one-dimensional array of x elements, all 
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initialized to zero. The ith element of the array contains the number of curvatures with magnitude i. For instance, 
the 24th element of the array for Pq is equal to 1, while the 25th element is equal to 0, since there are no curvatures 
equal to 25 in Pc- 

We use these arrays to generate the histograms in Section[3] Due to Matlab's memory constraints, we limit our 
Matlab arrays to 10 8 entries. So, to check for exceptions in the entire range of [10 6 ,5 • 10 8 ), we check each of the 
intervals [10 6 , 10 8 ), [10 8 ,2 ■ 10 8 ), . . . , [4- 10 8 ,5 • 10 8 ) individually. We have chosen to display the interval [10 6 , 10 8 ) 
in our figures in Section[3] 

To count primes less than x, we simply increment a sum whenever a prime curvature is produced. To count 
kissing primes less than x, we increment a sum whenever a prime curvature is produced and some other member 
of the curvature's quadruple is prime. 

It takes our algorithm 0(Np(x)) steps to compute Np(x), which is optimal since each node on the tree must be 
visited; that is, it is not possible to skip any quadruples. 

Our programs rely on Wayne and Sedgewicks's Stack data type and standard draw library ISWI . 



s 2^— (157, 397, 24, 28) 




Figure 10. The tree of quadruples for Pc, pictured up to 2 generations. 
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